Disentangling phonons from spins in ion-trap-based quantum spin simulators 
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We compute how phonon creation affects the fidelity of the quantum spin dynamics in trapped ion 
simulators. A rigorous treatment of the quantum dynamics is made by employing an exact operator 
factorization of the evolution operator. Although it is often assumed that phonon creation modifies 
the dynamics of the spin evolution, for an Ising spin-spin interaction in an external magnetic field, 
phonons have no effect on the probabilities of spin product states measured in the direction of 
the Ising model axis. Phonons play a much more important role in infiuencing the effective spin 
dynamics for Heisenberg or XY model spin simulators or for other observables, like witness operators 
in the Ising model. 
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In the early 1980's, Richard Feynman proposed that 
one could used controlled quantum-mechanical systems 
to simulate the many-body problem fl'|. In recent years, 
there has been significant success in trying to achieve 
this goal [2]. We focus here on one platform for per- 
forming analog quantum computation, the simulation of 
interacting quantum spins in a linear Paul trap with a 
finite number of ions [3|-|6|. In these systems, the clock 
states of the ions are the pseudospin states, which can 
be manipulated independently by a spin-dependent force 
driven by laser beams. The lasers couple the pseudospin 
states to the lattice vibrations of the trapped ions, which 
leads to effective spin-spin interactions when the phonon 
degrees of freedom are integrated out Recently, 
ion traps have been used to simulate the transverse-field 
Ising model with two 5] and three [3] ions in a trap. 
Experiments are usually run in the Lamb-Dicke regime 
where spontaneous phonon creation is low, but it is gen- 
erally believed that phonon creation, especially as one 
tunes the driving force close to resonance with one of 
the vibrational mode frequencies, will act to dephase and 
decohere the quantum dynamics leading to errors in the 
analog simulators that will build up as a function of time. 
Surprisingly, for the simplest class of emulators (namely 
the transverse-field Ising model) one can show that there 
is no dephasing or decoherence by the phonons to lowest 
order in the Lamb-Dicke parameter if one measures the 
probabilities of product states of spins where the spin 
quantization axis coincides with the Ising model axis. 
For more complicated models, like the XY or Heisen- 
berg models, phonon dephasing and decoherence cannot 
be removed, and hence this will create more stringent 
restrictions on the fidelities of quantum simulators for 
those models when simulated in ion traps. Moreover, 
when measuring entanglement witness operators, even in 
the transverse field Ising model, the phonons should af- 
fect the results. 

Quantum spin models are of interest for a number of 
reasons. Currently, there is a significant focus on the 
physics of frustration ^ , and how that can lead to the 



formation of a so-called spin liquid [lo| . These spin liq- 
uid states are complicated quantum many-body systems 
that exhibit significant entanglement of their wavefunc- 
tions, and could also exhibit emergent phenomena within 
their low-energy excitation spectra. Conventional com- 
putation fails to describe these systems well. Exact di- 
agonalization studies are limited to small size lattices, 
while quantum Monte Carlo simulations suffer from the 
sign problem and cannot reach low temperatures. Hence, 
they form an important test case for quantum simulators, 
and can allow the simulators to display interesting new 
physical phenomena. 

When N ions are placed in a (harmonic) linear Paul 
trap they form a nonuniform lattice, with increas- 
ing interparticle spacing as one moves from the cen- 
ter to the edge of the chain. The ions vibrate in all 
three spatial dimensions about these equilibrium posi- 
tions and are described by the phonon Hamiltonian 
Uph = Y^au^^au{aii,aai' + 1/2), in which a^^ is the 
phonon creation operator of the normal mode v along 
the spatial direction a ^ x,y, z. The ath spatial compo- 
nent of the jth ion's displacement operator J/Z" is related 
to the ath phonon normal mode amplitude (eigenvector 
of the dynamical matrix) 6"'' and the ath phonon oper- 
ator via SR^ = bf^Ku + al^yn/2ML0o,^ with M 
the mass of the ion. 

The two clock states used to describe the quantum 
spins are split by an energy difference Tiloq. A spin- 
dependent force on the ions is generated by using light to 
couple these _states to a third state via stimulated Raman 
transitions [isj . The light that couples to this state is de- 
tuned by a small frequency ^ called the detuning. The 
rotating wave approximation is used to map this three- 
state system onto a two-state system in the regime where 
ujo ^ fj, (these differ by at least 3 orders of magnitude 
in typical experiments). The two- level system is then 
represented by an effective laser- ion Hamiltonian isl - isj ] 
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where is the wavevector difference of the two inter- 
fering laser beams that generate the Raman coupUng. 
Here the Lamb-Dicke Umit Sk\S'Rj{t)\ ^ 1 is assmned so 
that the expansion for exp (— iJk ■ Slij) « 1 — iSk. ■ Stij 
is taken. The Roman index j denotes the different ions, 
flj is the Rabi frequency due to the Raman couphng at 
the jth ion site, and the spin operator trj is the Pauh 
spin matrix in the x-direction. While one can couple to 
any direction by tuning the phase difference of the Ra- 
man beams, for simplicity, and without loss of generality, 
we assume we have the lasers tuned to create a coupling 
to the spin component in the cc-direction only here. For 
example, we consider the ions crystallized along the z- 
axis of the trap, so the x- and y-directions correspond to 
transverse phonons. Then coupling the Raman lasers in 
the transverse direction minimizes effects of ion heating 
and allows for an identical spin axis for each ion Q. 

In experiments, one uses adiabatic quantum state evo- 
lution to evolve the ground state from an easily prepared 
initial state to the desired complex quantum state that 
will be studied. For spin models generated in an ion trap, 
it is easy to create a fully polarized ferromagnetic state 
via optical pumping, and then use a spin rotation to ori- 
ent the ferromagnetic state in any direction. Then, if one 
introduces a Hamiltonian with a magnetic field in the di- 
rection of the polarized state, it is in the ground state 
of the system. By slowly reducing the magnitude of the 
field and turning on the spin Hamiltonian of interest, we 
can reach the ground state of the Hamiltonian of interest 
(even in zero field). Hence, we need to add an addi- 
tional Zeeman term %B{t) = T^^=i ' with a time- 
dependent effective magnetic field B(i) (the coupling is 
to the different Pauli spin matrices). For example, a mag- 
netic field in the ^/-direction is made by driving a resonant 
radio-frequency field with frequency wq between the two 
hyperfine states to implement the spin flips [lij . The full 
Hamiltonian is then 'H{t) = Hph + 'HLi{t) +'HB{t)- 

Quantum dynamics for a time-dependent Hamiltonian 
are found by calculating the evolution operator as a time- 
ordered product U{t,to) — 7^exp[— i Jldt"Hit')/h] and 
operating it on the initial quantum state |V'(^o))- For 
the adiabatic evolution of the ground state, we start 
our system in a state with the spins aligned along the 
magnetic field and no phonons at time to- |'0(io)) — 
\0)ph ^ I tyty ■ ■ ■ ty)- In time-dependent perturbation 
theory, one rewrites the evolution operator in the inter- 
action picture with respect to the time-independent part 
of the Hamiltonian. This procedure produces a factorized 
evolution operator 



(2) 



which is the first step in our factorization procedure (the 
first factor on the left is called the phonon evolution 
operator). The second factor on the right is the evo- 
lution operator in the interaction picture which satis- 



fies an equation of motion given by i1ldUi{t,to)/dt = 
[Vi{t) + UBit)]Ui{tM). with Vi{t) = e^^\inph{t - 
to)/^'HLi{t)exp[-i'Hph{t ~ to)/n] since [hLB{t).,Uph\ = 
0. The only difference between 'H.Liii) and Vi{t) is 
that the phonon operators are replaced by their inter- 
action picture values: aa^, ^ aai, e:sq)[~iu!ai^{t — to)] and 
«L al^^eiip[iujai,{t ~ to)]- 

We now work to factorize the evolution operator fur- 
ther. Motivated by the classic problem on the driven 
harmonic oscillator in Landau and Lifshitz's Quantum 
Mechanics [lal and similar discussions in Gottfried's 



book [17[, we factorize the interaction picture evolution 
operator further Ui{t,to) = exp[—iWi{t)/h]U{t,to) with 
Wi{t) defined by Wi(t) = dt'Vi{t') (we call the factor 
on the right the phonon-spin evolution operator). One 
important result we use is that the multiple commutator 
satisfies [[Wi{t),Vi{t')],Vi{t")] = 0, which greatly sim- 
plifies the analysis below. 

The equation of motion for the spin evolution operator 
U{t, to) satisfies 
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in which the operator Hit) is given by the expression 
■H{t) = e^'^'^'^-ihdt+nB{t) + Vi{t)]e-^'^'^'l (4) 

The operator 'H(t) can then be expanded order by order 

as 

e^^^(*V/(Oe-^^^(*) = Vi{t) + ^[Wiit),Viit)] (5) 

N 

eiW'Ht)ij^(i)e-iM^Ht) = ^ |B(t) • a, (6) 

+ J[W,(0,B(i)- a,] + ...}, 

e^^^(*)*?i9te-^^^(*) = Viit) + ^^[Wiit),Vi{t)],{7) 

where we used the facts that dWi{t)/dt — Vi{t) and that 
[[Wj{t),Vj{t%Vj{tn]^0. 
Explicit calculations then yield 
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in which 7]^^ = Sax^k" ^/hj2MLU^ is the Lamb-Dicke 
parameter for the phonon mode with the ath compo- 
nent of the laser momentum Sk". 

We choose the magnetic field in the y-direction. Then 
the term j-[Wi{t'), HB{t')] is nonzero and proportional to 
spin operators in the z-direction times phonon operators. 
But it is also proportional to the Lamb-Dicke parameter 
Vaf ^ 1 which is small in current experiments. Hence, to 



lowest-order in the Lamb-Dicke regime, we drop all com- 
mutator terms with the magnetic- field piece of the Hamil- 
tonian (recall the magnetic field is a similar magnitude to 
the spin exchange parameters, so the additional Lamb- 
Dicke parameter multiplying the magnetic field makes 
it a less important term). The equation of motion can 
now be directly integrated. The spin evolution operator 
ij{t,to) becomes 



; dt' 
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which is the third factor for the evolution operator 
of the Ising model in a transverse field. The spin 
exchange terms Jjj' (t) — Jj^, + AJjj> (t) arise from 
the i[Wi (t) , Vi {t)] / 2ti commutator and include a time- 
independent exchange interaction between two ions Jj*^, 
and a time-dependent exchange interaction AJjji{t). 
The time-independent term determines the effective spin- 
spin Hamiltonian that is being simulated, while the time- 
dependent terms can be thought of as diabatic correc- 
tions, which are often small in current experimental set- 
ups, but need not be neglected. For simplicity, we set 
the initial time to = 0. Then the expression for the spin 
exchange interaction Jjj/ (t) is 

Jn'it) - 2^ yi j2 ' (12) 

X [lOxv — i^xu cos — 2/i sin uixi^i sin /ii] , 

which can be antiferromagnetic Jjj/ (t) > or ferromag- 
netic Jjji (t) < depending on the laser detuning /i and 



p^it) = E ••• E ■•• E m^phinxi,. 

n^i=0 n^„=0 n^H=0 



the detailed phonon mode properties w^i,, bj", and b^l^. 

We need to evaluate one more identity before 
arriving at our main result. We further fac- 
torize the entangled phonon-spin evolution operator 
exTp[-iWi{t)/h] into the product exp[-i r*j,(t)aj.^] 
exp[-« E,. ^xAt)axu] exp[- E^ Txuit)ri^it)/2] with the 
spin operator defined to be r*^(t) — Ej7f (*)'''j ^^'^ 
its complex conjugate is Txuit), in which the func- 
tion 7j''(t) satisfies jf'it) = njij^^b^'^/in^ - u^lS) x 
[exp{— ia;a;i/t}(jWa;i/ sin/Lt< -I- cos /it) — /i]. 

We are now ready to show our main result that 
phonons have no effect on the probabilities to observe 
any of the 2^ product states with a quantization axis 
along the Ising axis |/3) = | or ix) ®\'\xOt:\,x)®--- 
for the N ionic spins. Using the fundamental axiom of 
quantum mechanics, the probability Pp{t) to observe a 
product spin state |/3) starting initially from the phonon 
ground state (for all phonon modes) |0) and not measur- 
ing any of the final phonon states involves the trace over 
all possible final phonon configurations 



■ • ^xiy, ■ ■ ■ UxN \U{t,to)\0,0,...,0)ph<S>mO))\^ (13) 



where |$(0)) is any initial spin state (it need not be a 
product state). Note that since the pure phonon fac- 
tor of the evolution operator exp[—i7ipht/h] is a phase 
factor, it has no effect on the probabilities when eval- 
uated in the phonon number operator basis, so we can 
drop that factor. Next, the term exp[—i'^^Txi,(t)axv'] 
gives 1 when operating on the phonon vacuum to the 
right, so it can be dropped. We are thus left with 
three factors in the evolution operator. One involves 



exponentials of the phonon creation operator multiplied 
by spin operators, one involves products of spin opera- 
tors that resulted from the factorization of the coupled 
phonon-spin evolution operator factor, and one is the 
pure spin evolution factor U. The two factors that ap- 
pear on the left involve only Ising spin operators, and 
hence the product state basis is an eigenbasis for those 
operators. This fact allows us to directly evaluate the 
expression in Eq. ([T^. We expand the evolution of the 
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initial state at time t in terms of the product-state ba- 
sis m)) = IJitM = m%) = Y.p'Cp'{t)\fi') with 1/3') 
denoting each of the 2^ product state basis vectors and 
cpi [t) is a number. Using the fact that the product states 



satisfy the eigenvalue equation crj |/3') = to|^, |/3') with 
eigenvalues m|^, = +1 (for | fa;)) or —1 (for | ix))^ we 
arrive at the expression for the probability 



Pfi{t) = |c^(t)|'exp 
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We used the matrix element 
ph{nxi, ■ ■ ■ , nxvi ■ ■ ■ nxN\e 



''XL' 7 ■ 



-KM 



|0,...0)p^l5) 

-KNitM 



in the derivation. The summations become exponentials, 
which exactly cancel the exponential term and finally 
yield P/sit) = |c^(t)p, which is what we would have found 
if we evaluated the evolution of the spins using just the 
spin evolution operator U. Hence, phonons have no ob- 
servable effects on the probability of product states for the 
transverse-field quantum Ising model. If we do not mea- 
sure the probability of product states, then the terms 
from the coupled phonon-spin evolution operator remain 
spin operators, and one can show that the probabilities 
are changed by the phonons by simply evaluating the re- 
sult in Eq. (IT3l) directly. In other words, it is because the 
r operators are diagonal in the product space basis along 
the Ising axis that allows us to disentangle the phonon 
and spin dynamics. In cases where this cannot be done, 
we expect the phonon and spin dynamics to generically 
remain entangled. In particular, when one measures a 
typical entanglement witness operator, one measures the 
probabilities for the state to be projected onto product 
states that do not lie along the Ising axis. In this case, 
the phonons and spins will remain entangled, and the 
phonons will affect the expectation value of the witness 
operator. 

The simplest way to set up XY or Heisenberg model is 
to introduce additional laser beams with perpendicular 
i5k values so one gets pure spin-spin interactions along 
each of the different coordinate axes (if one cannot in- 
troduce perpendicular beams, then other methods can 
be tried, but they require strategies to reduce unwanted 
cross terms like a^a^). For the XY or Heisenberg model, 
we generate additional entangled spin-phonon evolution 
operator factors for the different coordinate directions 
{x and y for the XY model and all for the Heisenberg). 
But, since one cannot find a basis of quantum states that 
simultaneously diagonalizes the Pauli spin operators in 
more than one coordinate direction, the above derivation 
cannot be completed, and so the phonons will generically 
modify the probabilities for the spin states. 



In conclusion, we have explicitly shown that phonons 
do not affect the probability of observing spin prod- 
uct states in ion trap simulators of the transverse-field 
Ising model to lowest order in the Lamb-Dicke parame- 
ter. This implies that simulations of the transverse-field 
Ising model have the phonon motion and the spin evolu- 
tion disentangled, and hence will be much cleaner than 
simulations of the XY or Heisenberg model. Our results 
do not depend on the detuning lying far enough from res- 
onance, unless the proximity to resonance moves the sys- 
tem out of the Lamb-Dicke regime so higher-order terms 
need to be included, in which case phonon entanglement 
effects are likely to affect the results. 
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